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Using gauge/gravity duality, we study the creation and evolution of anisotropic, homogeneous 
strongly coupled A/" = 4 supersymmetric Yang-Mills plasma. In the dual gravitational description, 
this corresponds to horizon formation in a geometry driven to be anisotropic by a time-dependent 
change in boundary conditions. 



Introduction. — The realization that the quark-gluon 
plasma (QGP) produced at RHIC is strongly coupled 
[1] has prompted much interest in the study of strongly 
coupled non-Abelian plasmas. Hydrodynamic simula- 
tions of heavy ion collisions have demonstrated that the 
QGP produced at RHIC is well modeled by near-ideal 
hydrodynamics [5], which is a signature of a strongly 
coupled system. The success of hydrodynamic mod- 
eling of RHIC collisions suggests that the produced 
plasma locally isotropizes over a time scale Tjso ^ 1 
fm/c [5] . Understanding the dynamics responsible for 
such rapid isotropization in a far-from-equilibrium non- 
Abelian plasma is a challenge. 

Because of the difficulty in studying real time dynam- 
ics in QCD at strong coupling, it is useful to have a toy 
model where one can study the dynamics of a far from 
equilibrium, strongly coupled non-Abelian plasma in a 
controlled setting. One such toy model is strongly cou- 
pled A/" =4 supersymmetric Yang-Mills theory (SYM), 
where one can use gauge/gravity duality to study the the- 
ory in the limit of large and large 't Hooft coupling A. 
This has motivated much work devoted to studying var- 
ious non-equilibrium properties of thermal SYM plasma. 

We are interested in exploring the physics of isotropiza- 
tion in far-from-equilibrium non-Abelian plasmas, in the 
simplest setting which allows complete theoretical con- 
trol. This leads us to focus on the dynamics of ho- 
mogeneous, but anisotropic, states in strongly coupled, 
large Nc SYM. A conceptually simple way to create non- 
equilibrium states is to turn on time-dependent back- 
ground fields coupled to operators of interest. To cre- 
ate states in which the stress tensor is anisotropic, it is 
natural to consider the response of the theory to a time- 
dependent change in the spatial geometry. For simplic- 
ity, we limit attention to geometries which have spatial 
homogeneity (i.e., translation invariance in all spatial di- 
rections), an 0(2) rotation invariance, and a constant 
spatial volume element. The most general metric satis- 
fying these conditions may be written as 

d.s' = -dt" + e^"(*) dxl + e-2Bo(t) , (1) 

where x±^ = {xi,X2\- 

The function Bq (t) describes a time-dependent shear in 
the geometry; neglecting (4-dimensional) gravity, BQ{t) 



is a function one is free to choose arbitrarily. We will 
choose Boit) to be asymptotically constant as t — > ±oo. 
We will also choose the initial state to be the SYM vac- 
uum. A time-dependent geometry will do work on the 
quantum system. Consequently, the state in the distant 
future will be a non-vacuum state which (when the geom- 
etry is once again static) will be indistinguishable from 
a thermal state. During the evolution, because the met- 
ric ([ij changes in an anisotropic fashion, the resulting 
plasma will also be anisotropic with different pressures 
(i.e., stress tensor eigenvalues) in the longitudinal (a;||) 
and transverse {x±) directions. Spatial translation in- 
variance implies that no hydrodynamic modes can be ex- 
cited. Therefore, the non-equilibrium plasma produced 
by the changing metric ([iJ provides a nice laboratory 
to study the relaxation of non-hydrodynamic degrees of 
freedom in a far from equilibrium setting. We choose 

Bo(t) = ic[l -tanh(i/r)] . (2) 

For c ^ 0, this represents a time-dependent rescaling of 
lengths in transverse directions relative to those in the 
longitudinal direction, over a period of order r. The lack 
of any other scale in conformally invariant SYM implies 
that the final state energy density will be 0{t^'^). With- 
out loss of generality we measure all quantities in units 
where t — I. 

Gravitational description. — Gauge/gravity duality |3] 
provides a gravitational description of large SYM in 
which the 5d dual geometry is governed by Einstein's 
equations with a cosmological constant. Einstein's equa- 
tions imply that the boundary metric g^^{x), which char- 
acterizes the geometry of the spacetime boundary, is dy- 
namically unconstrained. The specification of the bound- 
ary metric serves as a boundary condition for the 5d Ein- 
stein equations. This reflects the fact that the dual field 
theory (which resides on the boundary) does not back- 
react on the boundary geometry, whereas the boundary 
geometry can influence the dynamics of the field theory. 

We consider a 5d geometry which coincides with AdSs 
in the distant past. This geometry is dual to the vacuum 
of SYM. A time dependent boundary metric g^u{x) will 
create gravitational radiation which propagates from the 
boundary into the bulk. This infalling gravitational radi- 
ation will lead to the formation of a horizon , which acts 
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as an absorber of gravitational radiation — any radiation 
which passes through the horizon cannot escape back to 
the boundary. At late times when the boundary geom- 
etry is no longer changing, the bulk geometry outside 
the horizon will relax and asymptotically become static. 
This is the gravitational description of thermalization in 
SYM. 

Diffeomorphism and translation invariance allows one 
to chose the metric ansatz 



dx\^ 



dx 



lU -r 2drdv , (3) 



where A, B, and E are all functions of the radial coor- 
dinate r and time v only. Infalling radial null geodesies 
have constant values of v (as well as x± and Out- 
going radial null geodesies satisfy dr/dv = \A. At the 
boundary, located at r — oo, the coordinate v coincides 
with the boundary time t. The geometry in the bulk at 
V > corresponds to the causal future of t > on the 
boundary. The form of the metric (|3| is invariant under 
the residual diffeomorphism r ^ r + f{v), where f{v) is 
an arbitrary function. 

With a metric of the form ^ , Einstein's equations may 
be reduced to the following set of differential equations: 






= E(S)' + 2E'S-2E^ 


(4a) 





= E(S)'-F ^{J:'B + B't) , 


(4b) 





= A" + 3B'B -12T.'±/T.^ +A, 


(4c) 





= ±+ 1{B^Y.~ A'±) , 


(4d) 





= E" + S , 


(4e) 



where 



for any function h(r,v), 
h' = drh, h = dyh - 



^Adrh. 



(5) 



Eqs. (4d) and (4el are constraint equations; the radial 



derivative of Eq. (4d) and the time derivative of Eq. (4e 



are implied by Eqs. (4a)-(4c 



The above set of differential equations must be solved 
subject to boundary conditions imposed at r = oo. The 
requisite condition is simply that the boundary metric 
g^^{x) coincide with our choice (jlj) of the id geometry. 
In particular, we must have 



lim S(r, v)/r — 1 , 



lim B{r,v) = Bo{v) . (6) 



One may fix the residual diffeomorphism invariance by 
demanding that 

lim \A{r, v) - r^l /r = . (7) 

r — *oo 

These boundary conditions, plus initial data satisfying 



the constraint (4el on some v — const, slice, uniquely 
specify the subsequent evolution of the geometry. 

Given a solution to Einstein's equations, the SYM 
stress tensor is determined by the near-boundary be- 
havior of the 5d metric . If denotes the gravi- 
tational action, then the SYM stress tensor is given by 

r^''(x) = (2/y=?JM) 6S0/S9UX) . 



Near the boundary one may solve Einstein's equations 
with a power series expansion in r. Specifically, A, B and 
E have asymptotic expansions of the form 



A{r,v) = ^ [&„(«) + a„(i;) log r]r2 " 

n=0 

B{r,v) = ^ [ 6„ (f ) -h/?„(i;) log r]r-", 

n=0 

S(r,u) = ^ [s„(w) -Ha„(i;) log r]ri-". 



(8a) 
(8b) 
(8c) 



n=0 



The boundary conditions ^ and Q imply that &o(w) = 
Bo{v), So(v) = 1, ao(w) = 1, and ai{v) = 0. Substitut- 
ing the above expansions into Einstein's equations and 
solving the resulting equations order by order in r, one 
finds that there is one undetermined coefficient, 64 (w). 
All other coefficients are determined by the boundary 
geometry, Einstein's equations, and &4(w) [10 . 

By substituting the above series expansions into the 
variation of the on-shell gravitational action, one may 
compute the expectation value of the stress tensor in 
terms of the expansion coefficients. This procedure has 
been carried out in Ref. [5j, so we simply quote the re- 
sults. In terms of the expansion coefficients, the SYM 
stress tensor reads 

Tt = {N'j27T')dmg{-£,V^,V^,Vii), (9) 
where (with b'^^'' = d^bo): 

21(&W)4 - 468(5[,'V^^o'^ 

(10b) 



hi 



1 

256 



hi - 



1 

768 



+ 10{b^Pr + 4bi'h^^^ + 6Abi'^ 
V\\ = - ia4-264 + 

+ 10(6(^^)2 + 4&W6(^)-1285{,^) 



21(6^)4 + 936(6«)25('^ 



(10c) 



Numerics. — One may solve the Einstein equations 
for the time derivatives t], B, and A". Define 



(4a|-(4c 



e(r,w) 
^{r,v) 



dw [S(w, v)'^ — hi{w, v)] — Hi{r, v) , 

(Ua) 

dw 



2e{w,v)B'{w,v)Y.{w,v)~^/^ 

-h2{w,v)\-H2ir,v), (lib) 
where is an indefinite (radial) integral of /i„. 



Then Eqs. (4al-(4c) are solved by 

± = -26 

B = 

A" = -4 - 246 E'E"-* + 



(12) 

(13a) 
(13b) 
(13c) 
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The functions hn{r,v) are not constrained by Einstein's 
equations — their presence inside the integrands of 



Eq. (Ill are compensated by the subtraction of their 
integrals i?„(r, w). However, since we have chosen the 



upper Hmit of integration in Eq. ( 11 1 to be r = cxd, the 



functions hn{r,v) must be suitably chosen so that the 
integrals (11) are convergent. The simplest choice to ac- 



complish this is to set hi{r,v) equal to the asymptotic 
expansion of Y.{r,v)^ up to order l/r*"', for some A: > 1, 
and to set equal to the asymptotic expansion of 

2Q{r,v)B'{r,v)/'S{r,v)^^^ up to order l/r*"'. In our nu- 
merical solutions reported below, we use fc > 4. This 
choice makes the large r contribution to the integrals in 



Eq. (Ill quite small. As the coefficients of the series ex- 
pansions ([8| only depend on &o(w) and 64(f) and their 
V derivatives, this choice determines hn{r,v) in terms of 
one unknown function 64 (w). 

With the subtraction functions hn specified by the 
aforementioned asymptotic expansions, integrating Eq. 



( 12 1 fixes the compensating integrals 77„ up to an in- 



tegration constant which is an arbitrary function of v. 



Integrating Eq. (13c I for A(r,v) introduces two further 



(v dependent) constants of integration. The most direct 
route for fixing these constants of integration is to match 



the large r behavior of the expressions (13a I and (13b) 
and the integrated version of Eq. ( 13c ) to the correspond- 



ing expressions obtained from the series expansions (|8|. 
This fixes all integration constants in terms of bo and 64. 

Our algorithm for solving the initial value problem 
with time dependent boundary conditions is as fol- 
lows. Given an initial geometry defined by _B(r, wq); 
one knows 64(^0). Integrating the constraint equation 
( [4e| , with the fully determined asymptotic behavior ( 8c I , 
yields S(r, wq). From this information, one can com- 

With A{r,vo), 



pute A(r, Wo) by integrating Eq. (13c 



B{r,V()) and S(r, ug) known, one can then compute the 
time derivative 9„i3(r, wo) from Eq. (13b) and step for- 
ward in time. 



B{r, Wo -I- Aw) w B{r, vq) + d^B{r, vq) Av . 



(14) 



Repeating the above process using this updated profile 
of B determines S and A at time vq + Av, from which 
one computes dyB for the next time step. For an initial 
geometry corresponding to the SYM vacuum, plus the 
choice ([2| of boundary data, one has 



B{r, 



T,{r,— 00) =r, A{r, —00) — r'^ . (15) 



An important practical matter is fixing the computa- 
tion domain in r — how far into the bulk does one want 
to compute the geometry? If a horizon forms, then one 
may excise the geometry inside the horizon as this re- 
gion is causally disconnected from the geometry outside 
the horizon. Furthermore, one must excise the geome- 
try to avoid singularities behind horizons [B] . To per- 
form the excision, one first identifies the location of an 



apparent horizon (an outermost marginally trapped sur- 
face) which, if it exists, must lie inside a true horizon 
[7] . We have chosen to make the incision slightly inside 
the location of the apparent horizon. For the metric (|3|, 
the location rh{v) of the ap paren t horizon is given by 
t{rh{v),v) = or, from Eq. |l3a| , e{rh{v),v) = . 

Results and Discussion. — Fig. [T] shows a plot of the 
energy density and transverse and longitudinal pressures 
produced by the changing boundary geometry ([!]), with 
c = 2. These quantities begin at zero in the distant past 
when the system is in its vacuum state, and at late times 
approach thermal equilibrium values given by 



T^i- ^{ir'N^TyS) diag(3, 1,1,1), 



(16) 



where T is the final equilibrium temperature. Non- 
monotonic behavior is seen when the boundary geometry 
changes most rapidly around time zero |llj . 
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FIG. 1: Energy density, longitudinal and transverse pressure, 
all divided by as a function of time for c = 2. 

Fig. |2] displays the congruence of outgoing radial null 
geodesies, for c ~ 2. The surface coloring shows A/r"^. 
In the SYM vacuum (i.e., at early times) this quantity 
equals 1, while at late times A/r'^ = 1 — (rh/r)'^- Ex- 
cised from the plot is a region of the geometry behind 
the apparent horizon. In the SYM vacuum, outgoing 
geodesies are given by l/r-t-w/2 = const., and appear as 
straight lines in the early part of Fig. |2]. In the vicin- 
ity of w = 0, when the boundary geometry is changing 
rapidly and producing infalling gravitational radiation, 
the geodesic congruence changes dramatically from the 
zero temperature form to a finite temperature form. As 
is evident from the figure, at late times some outgoing 
geodesies do escape to the boundary, while others fall 
into the bulk and never escape. Separating the 'escap- 
ing' and 'plunging' geodesies is one geodesic which does 
neither — this geodesic, shown as the black line in Fig.|2] 
defines the true event horizon of the geometry. 

Fig. [3] plots the area of the apparent and true event 
horizons, again for c — 2. Nearly all growth of the ap- 
parent horizon area occurs in the interval — 2 < w < 0, 
during which the boundary geometry is changing rapidly. 
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FIG. 2: The congruence of outgoing radial null geodesies. 
The surface coloring displays A/r'^. The excised region is 
beyond the apparent horizon, which is shown by the dashed 
green line. The geodesic shown as a solid black line is the 
event horizon; it separates geodesies which escape to the 
boundary from those which cannot escape. 
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FIG. 3: Area elements of the true event horizon and the 
apparent horizon as a function of time. 
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TABLE I: Final equilibrium temperature T and isotropization 
time Tiso (in units of or r), for various values of c. The 
isotropization time riao is the time at which the pressures 
deviate from their equilibrium values by less than 10%. 



time Tiso. (These quantities only depend on |c|.) We 
define Tiso as the time when the transverse and longi- 
tudinal pressures equal their final values to within 10%. 
When |c| > 2, we find that Tiso ~ 2t, while for |c| < 2, 
Tiso ~ 0.7/T. Since Tiso is only a few times longer than 
the time scale t over which the boundary geometry (jlj is 
changing, this measure of isotropization time should best 
be viewed as an upper bound on isotropization times as- 
sociated with plasma dynamics alone. Nevertheless, it 
is interesting to note that Tiso ~ 0.7/T corresponds to a 
time of ^ fm/c when T = 350 MeV, similar to estimates of 
thermalization times inferred from hydrodynamic mod- 
eling of RHIC collisions 17 . 
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In contrast, the area of the true horizon grows in the dis- 
tant past long before the boundary geometry is signifi- 
cantly perturbed. This is a reflection of the global nature 
of event horizons — the location of the event horizon de- 
pends on the entire history of the geometry. It has been 
argued |S] that it is the area element of the apparent 
horizon, pulled back to the boundary along v = const, 
infalling null geodesies, which should be identified with 
the entropy density (times 4Gjv) in the dual field theory. 

Table |T] shows, for various values of c, the final equilib- 
rium temperature T and a measure of the isotropization 



